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ABSTRACT 

Observations of the radio remnant of Supernova 1987 A using the Australia Telescope Compact Array in the 
12-mm band on 2003 July 31 (day 6002.7 after the explosion) give the first fully resolved radio image of the 
supernova remnant. The diffraction-limited image has a resolution of about 0."45, a factor of two better than 
that of the 3-cm images previously obtained. There is excellent agreement between the 12-mm image and a 
contemporaneous super-resolved 3-cm image. Super-resolution of the 12-mm image gives a further factor of 
two improvement in resolution, to 0"25, albeit with limited dynamic range. While the spatial distributions 
of the radio and X-ray emission are broadly similar, there are significant differences in detail with no corre- 
spondence in the regions of brightest emission. The 12-mm image is well modelled by a thick equatorial ring 
inclined at 43° to the line of sight. This, together with the common east- west asymmetry and the relatively 
steady increase in the radio flux density, suggests that the reverse shock is the main site for generation of the 
radio emission. 

Subject headings: supernova remnants — supernovae:individual(SN 1987A) — radio continuum: ISM 



1. INTRODUCTION 

Supernova 1987 A, the brightest and closest supernova (SN) 
in nearly 400 years, was first observed on 1987 February 24. 
The supernova, located in the Large Magellanic Cloud, illu- 
minated and ionized much of the circumstellar material, re- 
sulting in the dramatic optical rings seen in reco mbination- 
line emission (Burrows et alJl9 95llS ugerman et alJ2002) and 
scattered light (e.g., Crot ts et alJl995t lSugarman et al. 2005). 
Hubble Space Telescope imaging has revealed a series of 
"hot spots" around the inner edge of the central equatorial 
ring, indicating that the SN blast wave is beginning to im- 
pact on the denser circ umstellar material of the equatorial ring 
(Sugerman et al. 2002). Spectral imaging of high-velocity (~ 
12, 000 km s" 1 ) Lya and Ha emission has enabled mapping of 
the re verse shock formed behind the blast wave ( Michael et al. 
120031) . These observations show that in 2000 - 2001, around 
5000 days after the SN, the reverse shock had a radius of about 
75% that of the equatorial ring and was expanding at about 
80% of the blast-wave velocity. The shock appears about 40% 
brighter in the eastern quadrant and is concentrated within 30° 
of the equatorial plane. 

SN 1987 A was observed as a short-lived radio supernova in 
the weeks following the optical supernova ( Turtle et al. 1987). 
It was redetected in mid- 1990 (around day 1200) when in- 
creasing radio e mission signalled the birth o f the supernova 
remnant (SNR) JStavelev-Smith et alJ ll992). Observations 
with the Molonglo Observatory Synthesis Telescope and Aus- 
tralia Telescope Compact Array (ATCA) show that, at least up 
until the beginning of 2001, the SNR radio emission has in- 
creased nearly line arly with time at frequencies between 84 3 
MHz and 9 GHz <Ball et alJl200ll IManchester eTal]l2002l) . 
The radio emission appears to be optically thin synchrotron 
emission with a power-law spectrum of index about -0.9. 
From late- 1992 the SNR has been strong enough to resolve 
at 9 GHz, showing that radio emission is annular in form, 
contained within the optical equatorial ring and expanding at 
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about 3000 km s" 1 Jstavelev- Smith et alJl993llGaensler et alJ 
119971 Manches ter et alJl2002l) . The radio emission peaks on 
the eastern and western sides of the ring and is signficantly 
brighter on the eastern side. These images had an effec- 
tive resolution of about 0"5 and relied on the use of "super- 
resolution" techniques in which visibilities in the outer parts 
of the sampled u-v plane are given higher weight. The re- 
sulting images are inherently uncertain since the solution ex- 
trapolates into unsampled parts of the uv-plane. 

X-ray emission was first detecte d from SN 1987A about 
1500 days after the explosion ( Gorenstein et al J 119941 
iBeuermann et aljfl994l) . Recent observations w ith the Chan- 
dra X-ray Observatory JPark et alJl2004t 120051) show that the 
X-ray luminosity is increasing dramatically with a doubling 
time of about 600 da ys. Chandra im aging observations taken 
in 2000 December (Park et all 120021) showed different spa- 
tial distributions for the lower energy (0.3 — 1.2 keV) pho- 
tons and those at higher energies (1.2 - 8 keV). The lower 
energy emission was concentrated in hot spots around the in- 
ner edge of the equatorial ring and was closely correlated with 
the optical emission, whereas the higher energy emission was 
slightly more central and apparently related to the radio emis- 
sion with a pronounced east-west asymmetry. In the more re- 
cent images, this distinction appears to have disappeared with 
the emission more uniformly distributed around the ring at all 
energies. 

Receivers operating in the 12-mm band (16-26 GHz) were 
installed on three antennas of the ATCA in 2001 October and 
on the remaining three antennas in 2003 April, giving a 6-km 
maximum baseline. One of the major scientific justifications 
for development of the 12-mm system was to obtain radio im- 
ages of SNR 1987 A with a resolution comparable to those ob- 
tained by the Hubble Space Telescope and the Chandra X-ray 
Observatory . SNR 1987 A was observed at 12 mm on 2001 
October 27 (day 5360 after the SN) with a 1.5-km maximum 
baseline. While these observations just marginally resolved 
the SNR, they gave flux densit ies of 20 ± 1 and 18 ± 1 mJy a t 
17 and 19 GHz, respectively (IManchester & Gaensierfe OQl). 
In this paper we report on 12-mm observations of SNR 1987 A 
made on 2003 July 3 1 (day 6002.7) with six antennas over the 
full 6-km baseline in exceptionally stable atmospheric condi- 
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tions. 

2. OBSERVATIONS AND RESULTS 

SNR 1987A was observed for 1 1 hours on 2003 July 31 us- 
ing the 6D configuration of the ATCA with 12-mm receivers 
on each of the six antennas. Observations were obtained in 
two bands simultaneously, each of 128 MHz bandwidth, and 
centered on frequencies of 17.34 and 19.65 GHz respectively. 
Two phase calibrators were observed: 0541-7332 for 2 min 
at approximately 10-min intervals and 0516-621 for 2 min at 
approximately hourly intervals. Visibility phases for the sec- 
ondary calibrators were relatively stable on all baselines with 
maximum drift rates on 6-km baselines of about 50° per hour. 

All data processing was done using the MIRIAD 3 synthesis 
data an alysis system. Positions were referred to the Ma et al.l 
(1998) J2000 position for 0516-621: 05 h 16 m 44?9262(5), 
-62°07'05 ."398(5) (where the uncertainty in the last quoted 
digit is given in parentheses) derived from very long 
baseline interferometry (VLBI) observations. The de- 
rived J2000 position for 0541-7332 was 05 h 41 m 50?7669(8), 
-73°32'15."340(8). Flux densities were based on observa- 
tions of Mars which was assumed to have a disk temperature 
of 190 K, giving a flux density of 16.0 Jy at 17.34 GHz and 
20.5 Jy at 19.65 GHz on the date of the observations. Mars 
was resolved at all except very short baselines. The flux den- 
sity scale was transferred to the primary calibrator 0826-373 
using the MIRIAD task PLBOOT giving flux densities for this 
source of 1.75 Jy and 1.60 Jy at 17.34 and 19.65 GHz respec- 
tively. 

Figure^ shows the 12-mm image obtained by inverting 
the calibrated u — v data for both observed frequencies with 
uniform weighting. The resulting image was cleaned and 
restored with the diffraction-limited beam of half-maximum 
size 0"45 x 0"39 with the major axis at position angle 
2°(measured from north through east). Positions are absolute 
and th e cross marks t he position for SN 1987 A determined by 
iRevnolds et alJ (119951) by relating the VLBI and Hipparchos 
reference frames; this position has an uncertainty of about 
0"07. A super-resolved version of this image, obtained using 
the MAXEN routine and restored with a 0"25 circular beam 
is shown in Figure \T]p. Because of the relatively low signal- 
to-noise ratio of the 12-mm data, the dynamic range of this 
super-resolved image is only about 8:1. 

These observations give an integrated flux density for SNR 
1987A of 28.5 mJy and 24.5 mJy at 17.3 GHz and 19.6 GHz 
respectively, with an estimated uncertainty of 1 .5 mJy in each 
case. A 12-hour ATCA observation at 8.6 GHz on 2003 Au- 
gust 1 gave an integrated flux density for the SNR of 46 ± 1 
mJy. Taken at face value, these results imply a mean spectral 
index at day 6003 of -0.7 with an uncertainty of about 0.1. 

3. DISCUSSION 

Figure^ shows that the SNR ring structure is clearly re- 
solved with the diffraction-limited 1 2-mm imaging . Pre- 
viously, with the 3- cm imaging ( Gaensler et al.l H997t 
Manc hester et al.l2002l) . super-resolution techniques were re- 
quired to reveal this structure. While the consistency of suc- 
cessive super-resolved images of SNR 1987 A gave us confi- 
dence that they were reliable, a direct confirmation of that is 
valuable. This is provided by Figure |2] which shows an over- 
lay of the 12-mm diffraction-limited image on the 3-cm super- 
resolved image based on the ATCA observations of 2003 Au- 

3 See http://www.atnf.csiro.au/computing/software/miriad 



gust 1 . Although there is an excellent overall correlation be- 
tween the two images, there are slight differences. On both 
the eastern and western sides, the peak 3-cm intensities are 
displaced slightly southward of the 12-mm peaks. Whether or 
not this indicates different spectral indices in different parts of 
the ring remains to be confirmed by future observations. 

Figure |3] shows slices through the diffraction-limited 12- 
mm image (Figure^) at six position angles. These clearly 
show the pronounced east-west asymmetry with the eastern 
lobe stronger and centered at a slightly greater radial distance 
compared to the western lobe. It is also clear that the north- 
south diameter of the radio emission is less than the east-west 
diameter; based on the separation of the slice maxima, the 
ratio of diameters is about 0.75. Within the uncertainties, this 
ellipticity is consistent wit h the inclination of 43° derived for 
the optical equatorial ring (Crotts et al. 19951). 

Emission at the centre of the ring is significant and about 
20% of the average ring intensity. Images were generated 
within MIRIAD of two disks and two optically thin spheres of 
diameter 1 "55 and 1 "65 respectively and differenced to form 
a ring and shell of thickness 0" 1 in each case. These were then 
convolved with the diffraction-limited beam. Slices through 
these model images are shown in the upper part of Figure [3] 
Clearly, even a thin shell of emission has a larger surface 
brightness at its center than the observed SNR, whereas the 
thin ring or annulus has a lower (zero) surface brightness at 
its center. This suggests that something between a shell and 
an equator ial ring is the best m odel. The reverse shock de- 
scribed by Michael et al. (2003) has just this morphology. The 
top two curves in Figure|3]are slices at position angles 0° and 
90° through an optically-thin equatorial ring of latitude extent 
±30° and inner and outer radii of 0"6 and 0."8 respectively, 
inclined at 43° to the line of sight. Apart from the unmod- 
elled east-west asymmetry, these model slices agree remark- 
ably well with the corresponding cuts through the 12-mm im- 
age (lower part of Figure[3}. 

In previous analyses, the size and expansion rate of the ra- 
dio shell has been me asured by fitting an optically thin shell 
to the 3-cm u — v data (Gaensl er et alll997l: Manchest er et all 
120021) . The diameter obtained from such a fit to the 12-mm 
data is 1"649 ± 0."004, shown by the vertical dotted lines 
in Figure [3] The model slice for a thin shell shows that the 
emission maxima lie significantly inside the nominal shell ra- 
dius, consistent with the observed profiles. A similar fit to 
the 3-cm data of 2003 August 1 gives a shell diameter of 
1"652±0"003, fully consistent with the fit to the 12-mm 
data. 

Recent X-ray observations (Par klt^Lll200i r2005) show 
new hot spots emerging on the western side of the ring at both 
high and low energies, with the X-ray emission now form- 
ing a more complete ring. Figure |4] shows an overlay of the 
super-resolved 12-mm image aligne d by eye with a contem- 
poraneous broad-band X-ray image (Pa rk et al.l20 05). While 
the overall annular distribution of emission is very similar in 
the two images, in detail they differ substantially. The radio 
emission peaks on the eastern and western sides of the ring 
whereas the X-ray emission tends to peak in the "corners" of 
the ring with saddle points at the radio peaks. Hubble Space 
Telescope images of SNR 1987 A taken at a similar epoch in 
Ha light (reproduced in the review bv lMcCravll2005l) show 
a large number of hot spots distributed around the inner edge 
of the equatorial ring. Again, there is no detailed correspon- 
dence with either the radio or X-ray emission distribution. 

It is clear that the physical conditions under which the ra- 
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FIG. 1 . — Diffraction-limited (a) and super-resolved (b) images of SNR 1987 A at 12 mm wavelength from observations made on 2003 July 31 using the ATCA. 
The central cross marks the VLBI position of SN 1987 A and the flux density scale is in units of mJy/beam. For the diffraction-limited image, the off-source rms 
noise is approximately 120 /iJy. The half-maximum contour of the restoring beam is shown in the lower- left corner of each plot. 




RA offset (arcsec) 

FIG. 2. — Overlay of the diffraction-limited 12mm image (contours) on 
a super-resolved 3cm image from ATCA observations on 2003 August 1 
(greyscale). The contour intervals are -0.2, 0.2, 0.5, 1, 1.5, 2 and 3 mJy/beam 
and the greyscale wedge is labelled in units of mJy/beam. 



dio emission, optical (Ha) hot spots and the X-ray emission 
are generated are different, although all result from the inter- 
action of the expanding supernova ejecta with circumstellar 
material. The Ha hot spots are clearly associated with the in- 
ner edge of the equatorial ring and have little or no azimuthal 
asymmetry. The radio emission is concentrated on the east- 
ern and western sides and has a pronounced east-west asym- 
metry, with a much brighter eastern lobe. Initially the X-ray 
emission shared some of these features but, with the recent 
rapid brightening, it more resembles the optical distribution. 
The time evolution of the intensity of the optical and X-ray 
emission is also similar, being very non-linear with a rapidly 
increasing gradient, whereas the radio emission has increased 
almost linearly over the past decade with j ust a very small 

increase in gradient (iManchester et al.l2002l) . 

As discussed by previous author s (e.g., iSugerman et al.l 
l2002HPark et al.l2004l:lMcCravl2005L and references therein), 
these results suggest that the optical emission and at least the 
rapidly increasing X-ray emission are generated by the inter- 
action of the supernova blast wave with dense regions of the 
equatorial ring which extend inward. The fact that the emis- 
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FIG. 3. — Slices through the 12-mm diffraction-limited image at six posi- 
tion angles. Dashed lines show the zero for each slice. The offset is radial 
from the VLBI position and is positive toward north and/or east. The verti- 
cal dotted lines show the radius of a thin shell fitted to the u - V data. Slices 
through a thin ring, a thin shell and an inclined thick ring are shown in the 
upper part of the figure - see text for details. 



sion is now increasing rapidly and becoming more uniform 
around the ring suggests that the blast wave is now encounter- 
ing the main ring structure. Furthermore, the similar spatial 
distribution and similarly short evolution timescales shows 
that the optical and X-ray emitting regions coexist behind the 
advancing shock front and are generated contemporaneously. 

In contrast, the radio emission began to increase at an ear- 
lier time and has increased more-or-less steadily since. This 
suggests that the relativistic electrons responsible for the syn- 
chrotron emission are more widely distributed behind the 
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FIG. 4. — Overlay of the super-resolved 12mm image (contours) on a Chan- 
dra 0.3 - 8 keV image taken on 2003 July 8 (greyscale). The contour intervals 
are -0.2, 0.2, 0.5, 1, 1.5 and 2 mJy/beam and the greyscale is linear in X-ray 
counts after deconvolution and smoothing. 



shock front, perhaps filling the region between the forward 
and reverse shocks. It is possible that the reverse shock is the 
main site for acceleration of the synchrotron-emitting elec- 
trons. This idea is supported by the east-west asymmetry 
common to the radio emission and the reverse shock and by 
excellent agreement of the model slices in Figure [3]based on 
this morphology with the corresponding slices through the ob- 
served image. The lifetime of the sychrotron electrons is very 
long. Hence, unless adiabatic expansion or other loss pro- 
cesses are important, the rate of generation of synchrotron- 
emitting electrons has been relatively stable since day 1200 
when the radio emission was first detected. At least at recent 
times, the area of the shock surface (whether it is the forward 



or the reverse shock) has been growing roughly as f 2 , sug- 
gesting that the rate of particle acceleration per unit area of 
the shock front has been declining at roughly this rate. 

The measured flux densities of the SNR at 17 and 19 GHz 
are above the values extrapolated from lower frequencies. 
The SNR spectrum measured between 0.8 and 10 GHz is be- 
coming flatter, with a g radient of roughly 4 x 10~ 5 per day 
jManchester et alJ|2002|). The observed slope at day 6000 is 
about -0.8 (Stavelev-Smith et al. 2005). This is marginally 
consistent with the value derived from the measured 12-mm 
and 3-cm flux densities at this time, -0.7 ±0.1, suggesting 
that the rate of spectral flattening may be increasing or that 
the radio spectrum may not be power-law over the whole 
range from 0.8 to 19 GHz. We note that diffusive shock ac- 
celeration theories predict a concave synchrotron spectrum 
and there is some evidence for this in o ther SNRs (e.g., 
iRevnolds & Ellison! 1992tlJones et all2003l) . 

The increasing radio flux density of SNR 1987 A will soon 
give super-resolved images with higher resolution and im- 
proved dynamic range, allowing a more detailed comparison 
with the optical and X-ray images. This should elucidate the 
relationship between the emission processes in the various 
bands, showing whether or not the acceleration of the radio 
sychrotron electrons occurs at the reverse shock as suggested 
above. The increasing flux density will also allow a more pre- 
cise measurement of the frequency and spatial dependence of 
the radio spectral index. 
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